THE NON-SINGLET QCD EVOLUTION KERNELS IMPROVED BY 
RENORMALON CHAIN CONTRIBUTIONS 



S. V. Mikhailovf] 

Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical 

Physics, 

14 1980, Moscow Region, Dubna, Russia 

'66 ' 

^ ■ Abstract 

Os ; 

Closed expressions are presented for the contributions to QCD non-singlet forward evolu- 
tion kernels P{z) for the DGLAP equation and to V{x, y) for non-forward (ER-BL) evolution 
I equation for a certain class of diagrams which include renormalon chains. Calculations are 

performed in covariant ^-gauge, in the MS scheme. The assumption of "naive nonabelian- 
ization" approximation for kernel calculations is discussed, and a special choice of the gauge 
parameter ^ = — 3 is analyzed in this context. Partial solutions to the ER-BL evolution 
equation are obtained. 
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1 Introduction 



Evolution kernels are the main ingredients of the well-known evolution equations for parton distri- 
bution of DIS processes and for parton wave functions in hard exclusive reactions. These equations 
describe the dependence of parton distribution functions and parton wave functions on the renor- 
malization parameter /x^. Here I discuss the diagrammatic analysis and multiloop calculation of 
the forward DGLAP evolution kernel P{z) and non-forward Efremov-Radyushkin-Brodsky- 
Lepage (ER-BL) kernel V{x,y) in a class of "all-order" approximation of the perturbative 
QCD. The regular method of calculation and resummation of certain classes of diagrams for these 
kernels has been suggested in 0. These diagrams include the chains of one- loop self-energy 
parts (renormalon chains) into the one-loop diagrams (see Fig. 1). Here the results for both the 
kinds of kernels (DGLAP and ER-BL), obtained earlier in the framework of a scalar model in 

six dimensions with the Lagrangian Lint = g^^{4'ii^iV)(Q) with the scalar "quark" flavours (ipi) 

i 

and "gluon" {(p), are extended to the non-singlet QCD kernels. For the readers convenience some 
important results of the paper would be reminded. 

The insertion of the chain into "gluon" line ("chain-1" in ^) of the diagram in Fig.l a,b and 
resummation over all bubbles lead to the transformation of the one-loop kernel ( see, e.g., |Q) 
aPo(-2) = az = a(l — z) into the kernel P^^\z; A) 



T-, / \ — chain— 1 

ar'Q[z) = az — > 



P^^\z;A) = az 



izr^ii 



A) 



7^(0) 



where A = aNf'j^{0), a 



(4vr) 
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Here, 7^(^)(£) are the one-loop coefficients of the anomalous dimensions of quark (gluon at Nf = 1) 
fields in D-dimension {D = 6—2e) discussed in 0; for the scalar model 7v>(£^) = 7i/p(^) = B{2—e, 2 — 
e)C{e), and C{e) is a scheme-dependent factor corresponding to a certain choice of an MS -like 
scheme. The argument A of the function '^^p{A) in (|l|) is the standard anomalous dimension (AD) 
of a gluon field. So, one can conclude that the "all-order" result in (|lD is completely determined 
by the single quark bubble diagram. The resummation of this "chain-1" subseries into an analytic 
function in A shouldn't be taken by surprise. Really, the considered problem can be connected with 
the calculation of large Nf asymptotics of the AD's in order of '\-/Nf. An approach was suggested 
by A. Vasil'ev and collaborators at the beginning of 80' to calculate the renormalization-group 
functions in this limit, they used the conformal properties of the theory at the critical point g = Qc 
corresponding to the non-trivial zero Qc of the D-dimensional /3-function. This approach has been 
extended by J. Gracey for calculation of the AD's of the composite operators of DIS in QCD in any 
order n of PT, 0. I have used another approach, which is close to 0; contrary to the large Nf 
asymptotic method, it does not appeal to the value of parameters NfTji, Ca/'^ or Cp, associated 
in QCD with different kinds of loops. To illustrate this feature, let us consider the insertions of 
chains of one-loop self-energy parts into the "quark" line of diagram Fig. la ("chain-2" in 0). 
Contributions of these diagrams, calculated in the framework of the above scalar model, do not 
contain the parameter Nf and can be summarized into the kernel P^'^\z; B) |Q 



aPo{z) 



az 



chain— 2 



P^^\z-B) 



az 



1+5- 



dB 



where B = 07^(0), 



according to the same approach. Following this way, the "improved" QCD kernel P^^\z; A) has 
been obtained in for the case of quark or gluon bubble chain insertions in the Feynman gauge 
and in M for general case of the mixed insertions in ^- gauge. 
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In this talk, we present the QCD results similar to Eq.(|I]), for each type of diagrams appearing 
in the covariant ^- gauge for the DGLAP non-singlet kernel P{z; A). The analytic properties of the 
function P{z; A) in variable A are analyzed. The assumption of "Naive Nonabelianization" (NNA) 
approximation p for the kernel calculation [1^ is discussed and its deficiency is demonstrated. 
The ER-BL evolution kernel V{x, y) is obtained in the same approximation as the DGLAP kernel, 
by using the exact relations between P and V kernels |Tl|, |^ for a class of "triangular diagrams" 
in Fig. 

m 



1. The considered class of diagrams represents the leading Nf contributions to both kinds 
of kernels. The partial solutions for the ER-BL equation are derived (compare with 
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2 Triangular diagrams for the DGLAP evolution kernel 

Here, the results of the bubble chain resummation for QCD diagrams in Fig.l a,b,c for the DGLAP 
kernel are discussed. These classes of diagrams generate contributions ~ (a^ ln[l/2;])" in any 
order n of PT. Based on the resummation method of Ref. in the QCD version, one can derive 
the kernels P^^"-''^'^^ (corresponding to the diagrams in Fig.l) in the covariant ^— gaugeQ 



6{x — kn) 




MC 




MC 




(b) 



(d) 



+ - 



+ 



+ 



Figure 1: The diagrams in figs, la - Ic are the "triangular" diagrams for the QCD DGLAP 
kernel; dashed line for gluons, solid line for quarks; black circle denotes the sum of all kinds of the 
one-loop insertions (dashed circle), both quark and gluon (ghost) or mixed chains; MC denotes 
the mirror-conjugate diagram; Id is an example of diagram for the non-forward ER-BL kernel. 



P^^''\z;A) 
P'^^^\z;A) 
P^^''\z- A) 



a.Cp'lz ■ (1 - Afz 



I 7,(0) 



ttsCp ■ S{1 



7,(0) 



,{l-A)^,{A) 



2z 



l-A 



ttsCp ■ 5{1 — z 



i-^7,(^);+' 

A(3-2A) 7g(0) 



,(2-yl)(l-A)7,(A) 



(3) 
(4) 
(5) 



where as 
and A = 



as 
4^' 



-a 



Cf = {Nl - l)/2Nc, Ca = and Tr = - are the Casimirs of SU(Ar,) group, 
i7g(0). The function 7g(£:) is the one-loop coefficient of the anomalous dimension of 



^The gauge parameter ^ is defined via the gluon propagator in lowest order iD^i^{k^ 



9^.u + (e - 1) 
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gluon field in D-dimension, here D = 4: — 2e. In other words, it is the coefficient Zi{e) of a simple 
pole in the expansion of the gluon field renormalization constant Z, that includes both its finite 
part and all the powers of the ^-expansion. The function 7g(e) thus defined is an analytic function 
in the variable e by construction. Equations (0) - are valid for any kind of insertions, i.e., 
jg = 7^'^^ for the quark loop, jg = 7^^^ for the gluon (ghost) loop, or for their sum 

^g{A,0 = i'\A) + ^^f\A,0; 

when both kinds of insertions are taken into account. The sum of contributions (^, (^), 
in P^^\z; A,^) which has the expected "plus form" 



zz 



ttsPoiz) 



+ 



-Af 

2z 
1 - z 



2z 



l-A 



results 

(6) 
(7) 



here, for comparison with (^, the one- loop result asPo{z) is written down, the latter can be 
obtained as the limit P^^^ (2;; A — »• 0, 0- Note that in (P) the 5{l—z) - terms are exactly accumulated 
in the form of the [...]+ prescription, and the ^ - terms successfully cancel. This is due to the 
evident current conservation for the case of quark bubble insertions, including the gluon bubbles 
into consideration merely modifies the effective AD, 7g(y4,^) 'y^^\A), conserving the structure 



of the result (y) , see 0, 
gluon (ghost) loops (see. 



|]. Substituting the well-known expressions for 7g(e) from the quark or 
e.g., m) 

= -mfTRB{Dl2,D/2)C{e) 



.0 



B{D/2 - l,D/2 




(9) 



into the general formulae (^ - (^, and (^ one can obtain P^^\z]A,^) for both the quark and 
the gluon loop insertions simultaneously. Here, the coefficient C{e) = r(l — e)T{l + e) implies 
a certain choice of the MS scheme where every loop integral is multiplied by the scheme factor 
T{D/2 — 1) {fj!^ / Any . The renormalization scheme dependence of P^^\z; A) is accumulated by the 
factor C{e) |. Of course, the final result (P) will be gauge- dependent in virtue of the evident gauge 



dependence of the gluon loop contribution 7^^)(£, ^), in this case, e.g. 



A{0 = -aslMO 




Ca - -NfTn 



(10) 



is the contribution to the one-loop renormalization of the gluon field. The positions of zeros of 
the function 7g(A, ^) in A, which manifest the poles of P{z;A,E,), also depend on ^. The kernel 
P'^^\z; A,^) became gauge-invariant in the case when only the quark insertions are involved, i.e., 

7^^) 7,; A = A^'') = -as^f\Q) = as^TnNf, and P^'\z; A^'^^) ^ P^'\z;A,0 as it was 

_ o 
presented in [y]. It is instructive to consider this case in detail. To this end, let us choose the 

common factor 7^'^'*(0)/7^'^''(A) in formula (|^) for the crude measure of modification of the kernel in 

comparison with the one-loop result asPo{z). Considering the curve of this factor in the argument 

A in Fig.2, one may conclude: 

■^For another popular definition of a minimal scheme, when a scheme factor is chosen as exp(c • e) , c = —75 + . . . 
instead of r{D/2 — 1), the coefficient C{e) does not contain any scheme "traces" in final expressions for the 
renormalization-group functions. 
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Figure 2: The curve of the factor 7^(0) /7c, (A), the arrow on the picture corresponds to the point 
A = l/vr {i.e. = 6/Nf), the first singularity appears at A = Aq = 5/2. 



(A) 



(i) the range of convergence of PT series corresponds to the left zero of the function 7^ 
and is equal to Aq = 5/2, that corresponds to a° = 15ii/Nf, so, this range looks very broad [J, 
as < 571 at Nf = 3; 

(ii) in spite of a wide range of PT fidelity, the resummation into Pjj^^ (z; A) is substantial - 
two zeros of the function P^^^ (z; A) in A appear within the range of convergence (it depends on a 
certain MS scheme); 

(iii) the factor Jg'^^O) / 'Jg'^\A) decays quickly with the growth of the argument A. Really, if 
we take the naive boundary of the standard PT applicability, = 1 at A'^^ = 3, A^'^^ = l/(27r), 
then this factor falls approximately to 0.7, at Nf = 6, A^'^^ = l/vr it falls to 0.5, see arrow in Fig. 
2; thus, the resummation is numerically important in this range. 

Note at the end that Eq.(|^) could not provide valid asymptotic behavior of the kernels for 
z 0. A similar 2;-behavior is determined by the double- logarithmic corrections which are most 
singular at zero, like (^a5ln^[z]^ \T^]. These contributions appear due to renormalization of 
the composite operator in the diagrams by ladder graphs, etc. rather than by the triangular ones. 



3 Analysis of the NNA assumption for kernel calculations 

The expansion of Pj^^^z; A) in A provides the leading {asNf ln[l/2;])'^ dependence of the kernels 
with a large number Nf in any order n of PT But these contributions do not numerically 
dominate for real numbers of flavours Nf = 4, 5, 6. That may be verified by comparing the 
total numerical results for the 2- and 3-loop AD's of composite operators (ADCO) presented 
T6[| with their A^/-leading terms (see ADCO in Table 1). Therefore, to obtain a satisfactory 



m 



agreement at least with the second order results, one should take into account the contribution 
from subleading A^/-terms. As a first step, let us consider the contribution from the completed 
renormalization of the gluon line - it should generate a part of subleading terms. Below we 



'*Here we consider the evolution kernel P{z, A) by itself. We take out of the scope that the factorization scale /i^ 
of hard processes would be chosen large enough, /i^ > to^, where the p-meson mass nip represents the characteristic 
hadronic scale. Following this reason, the used coupling as{fi'^) could not be too large. 
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shall examine an exceptional choice of the gauge parameter ^ = —3. For this gauge the coefficient 
of one-loop gluon AD 73(0, —3) coincides with bo, the one-loop coefficient of the /3-function f\ 
Therefore this gauge may be used for a reformulation of the so-called NNA proposition to 



kernel calculations. Note, just this value of ^ has been used in [|T5| to estimate the total gluon 
contribution only from the gluon bubble in order of to the process of e"*" e~ annihilation. 

Table 1. 

The results of r(i 2)(^) calculations ( r(n) = dxx'^P{z)) performed in different ways, exact nu- 
merical results from [15] and approximation obtained from P{z, A,C,) with C, = —3; both numerical 
and analytical exact results are emphasized by the bold print. 





r(i)(n) 


r(2)H 




CfCa 


NfCp 




Nf ■ CfCa 




n=2 












Exact 


13.9 


-2.3704 


86.1 + 21.3 C(3) 


-12.9 - 21.3 C(3) 


-0.9218 


e = -3 


11.3 




-42.0 


12.9 




n=4: 












Exact 


23.9 


-4.9152 


140.0 + 19.2 C(3) 


-18.1 - 41.9 C(3) 


-1.5814 


e = -3 


23.5 




-76.0 


23. 




n=6 












Exact 


29.7 


-6.4719 


173+ 19.01 C(3) 


-20.4 - 54.0 C(3) 


-1.9279 


e = -3 


31.1 




-95.6 


28.5 




n=8 
Exact 

e=-3 


33.9 

36.3 


-7.6094 


196.9 + 18.98 C(3) 

-109.0 


-21.9 - 62.7 C(3) 

32.3 


-2.1619 


n=10 












Exact 

e = -3 


37.27 

41.00 


-8.5095 


216.0+ 18.96 C(3) 

-119.28 


-23.2 - 69.6 C(3) 

35.24 


-2.3366 


n=12 












Exact 

e=-3 


40.02 

44.64 


-9.2555 


7 

-127.61 


7 

37.58 


-2.4753 



the 



To obtain the NNA result in a usual way, one should substitute the coefficient 6o for Ig"^^ (0) into 
expression for A^'^'' by hand (see, e.g., fl^). Note, the use of such an NNA procedure to improve 



Here, for the /3(as)-function we adapt /3(as) — —boo^ 



bo = ^Ca - ^NfTR 
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pW(z;A) leads to poor results even for Pi{z) term of the expansion; a similar observation was 



also done in |jT^ . The NNA trick expresses common hope that the main logarithmic contribution 
may follow from the renormalization of the coupling constant. This renormalization appears as a 
sum of contributions from all the sources of renormalization of , corresponding diagrammatic 



analysis for two-loop kernels is presented in ||TT| , ^. In the case of the ^ = —3 gauge the one- loop 



gluon renormalization imitates the contributions from these other sources and the coefficient 60 
appears naturally. 

The expansion of kernel P^^\z; A, ^) generates partial kernels aj,P(^i-^(z), alP(2){z), . . . which in 
their turn produce ADCO r(i)(n), r(2)(n), . . . according to the relation T{n) = Jq dzz'^P{z). 



These elements of ADCO and a few numerical exact results from |16] are collected in Table 1, let 



us compare them: 

(i) we consider there the contribution to the coefficient T(i-){n) which is generated by the gluon 
loops and associated with Casimirs CfCa/'^i the C|^-term is missed, but its contribution is 
numerically insignificant. It is seen that in this order the C7i?CA-terms are rather close to exact 
values (the accuracy is about 10% for n > 2) and our approximation works rather well; 

(ii) in the next order the contributions to r(2)(n) associated with the coefficients Nf ■ CfCa 
and C^Cp are generated, while the terms with the coefficients Cp, Nf ■ Cp, CpCA are missed. In 
the third order, contrary to the previous item, all the generated terms are opposite in sign to the 
exact values, and the = —3 approximation" doesn't work at all. 

So, we need the next step to improve the agreement - to obtain the subleading A^/-terms by 
the exact calculation. 



4 The non-forward ER-BL evolution kernel 



Here we present the results of the bubble resummation for the ER-BL kernel V{x,y). It can be 
derived in the same manner as it was done for the DGLAP kernel P{z), see Appendix A in ||^. 
On the other hand V{x, y) can be obtained as a "byproduct" of the previous results for P{z), i.e., 
we shall use again P] the exact relations between the V and P kernels established in any order of 



PT |TT[] for triangular diagrams. These relations were obtained by comparing counterterms for the 
same triangular diagrams considered in "forward", Fig.la, and "nonforward" , Fig. Id, kinematics. 

Let the diagram in Fig.la have a contribution to the DGLAP kernel in the form P{z) = 
p{z) + 6{1 — z) ■ C; then its contribution to the ER-BL kernel (Fig. Id) is 



V{x,y)=C[e{y>x) 



'—dz + 6{y — x) ■ C, 



:il) 



where C = 1 + (x ^ x, y — ?/) to take into account the mirror-conjugate diagram. From relation 
([TT|) and Eqs. (|^), (^ for p(^"'^) we immediately derive the expression for the sum of contributions 

Y(la+lc) 



V^^''+^''\x,y-A,i) = asCp2-C 



7g(0,O 



(12) 



that may naturally be represented in the "plus form". Expression (|T2D can be independently 
verified by other relations reducing any V to P [Tl|, |18| (see formulae for the V ^ P reduction 
there) and we came back to the same Eqs.(|^), (^ for p(^"-'^\ Moreover, the first terms of the 
Taylor expansion of V^^°'''^\x,y; A) in A coincide with the results of the two-loop calculation in 
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||rT|. The relation P ^ V similar to Eq.(|TT]) has also been derived for the diagram in Fig. lb 



therefore, substituting Eq.(^ into (|l^) we obtain 



e(y > x)^P^^''^ 



X 



V'^^'\x,y-A,i)=a,CF2-C 



l-A 



e{y > x) ( - 



y -X 



(13) 



(14) 



Collecting the results in (|T2D and (|T^ we arrive at the final expression for V^^^ in the "main 
bubbles" approximation 

l-A / 



V^'\x,y;A,0 = asCF2-C 



X 



9{y >x)[- 

\y 



A + 



y-x^ 



(15) 



which has a "plus form" again due to the vector current conservation. The contribution V^^^ in 
(|T5|) should dominate for Nf ^ 1 in the kernel V. Besides, the function V^^\x, y; A,^) possesses an 
important symmetry of its arguments x and y. Indeed, the function V{x, y;A,^) = V^-^^ {x, y; A,^)- 
[yyY^^ is symmetrical under the change x ^ y, V{x,y) = V{y,x). This symmetry allows us to 
obtain the eigenfunctions ipni^) of the "reduced" evolution equation Q 
1 

V^'\x, y; A)My; A)dy = T{n; A)^„(x; A), (16) 



My, A) ~ (yy)'^^''^'"^ Cy^'Ky - y), here d^A) = {Da - l)/2, = 4 - 2A, (17) 

and d^{A) is the effective dimension of the quark field when the AD A is taken into account; 
C^\z) are the Gegenbauer polynomials of an order of a. The partial solutions $(x; as, I) of the 
original ER-BL-equation ( where / = ln{ fi"^ / fi^)) 

{fi%2 + P{as)da,) $(x; as, I) = V^^\x, y; A) $(y; a„ l)dy (18) 

are proportional to these eigenfunctions 'ipn{,x; A) for the special case [3{as) = 0, see, e.g. 

In the general case /3(as) ^ let us start with an ansatz for the partial solution of Eq. (p!8D , 
^n{x;as,l) ~ Xn{cLs,l) -ipnix-jA), and the boundary condition is Xn{O's,0) = 1; $„(a;;as,0) ~ 
ipn{x; A). For this ansatz Eq.(p!8|) reduces to 

{fi^d^2 + (3{as)da,) In (<l>„(x; a„ /)) = T{n; A). (19) 

In the case n = the AD of the vector current r(0; A) = 0, and the solution of the homogeneous 
equation in (|T^) provides the "asymptotic wave function" 

%ix; as, I) = A) ~ ((1 - x)xY'-^\ (20) 

where A = —as{fi'^)'y{0,C,) and aslfi"^) is the running coupling corresponding to (3{as). A similar 
solution has been discussed in in the framework of the standard NNA approximation. Solving 
simultaneously Eq. (|1^) and the renormalization-group equation for the coupling constant ds we 
arrive at the partial solution $„(a;; ds, I) in the form 

f^sit^^) T{n,A) 

/ 

Ja^f^l) /3(a) 

Recently, a form of the solution ~ ipn{x; A) with A = —asbo has been confirmed in [|12| by the 
consideration of conformal constraints [|19| on the meson wave functions in the limit Nf ^ 1. 



$„(x, as) ~ Xn(/i^) • i'n{x; A); where Xn(/i^ 



exp 



-da 



(21) 
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5 Conclusion 



In this paper, I present closed expressions in the "all order" approximation for the DGLAP kernel 
P{z) and ER-BL kernel V{x, y) appearing as a result of the resummation of a certain class of 
QCD diagrams with the renormalon chain insertions. The contributions from these diagrams, 
P^^\z;A) and V^^\z; A), give the leading Nj dependence of the kernels for a large number of 
flavours Nf ^ 1. These "improved" kernels are generating functions to obtain contributions to 
partial kernels like a^""*'^''P(„)(2;) in any order n of perturbation expansion. Here A ~ is a new 
expansion parameter that coincides (in magnitude) with the anomalous dimension of the gluon 
field. On the other hand, the method of calculation suggested in does not depend on the 
nature of self-energy insertions and does not appeal to the value of the parameters NfT^, Ca/'^ 
or Cp associated with different loops. This allows us to obtain contributions from chains with 
different kinds of self-energy insertions, both quark and gluon (ghost) loops, see The prize 
for this generalization is gauge dependence of the final results for P^^\z] A) and V^^\z] A) on the 
gauge parameter ^. 

The result for the DGLAP non-singlet kernel P'^^^z; A(^), ^) is presented in (^ in the covariant 
^-gauge. The analytic properties of this kernel in the variable are discussed for quark bubble 
chains only, and in the general case for an exceptional gauge parameter = —3. For the latter 
case, P^^\z] A{—3), —3) reproduces two-loop anomalous dimensions a^r(i)(n) with good accuracy, 
while the standard "naive nonabelianization" proposition fails at this level. But in the third order 
in tts the "^ = —3 approximation" is insufficient, see quantities r(2)(r;,) in Table 1. 

The contribution V^^\x,y; A{^),^) to the non-forward ER-BL kernel (p!5|) is obtained for the 
same classes of diagrams as a "byproduct" of the previous technique [Bl [Tl| . The partial solutions 



(|20| , ^Tj) to the ER-BL equation are derived. 

The obtained results are certainly useful for an independent check of complicated computer 
calculations in higher orders of perturbation theory (PT), similar to |]TB[; they may be a starting 
point for further approximation procedures. 
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